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1. INTRODUCTION 

The power quality (PQ) term originates from the idea that the AC voltage supplied to the domestic 
and industrial equipment must be unpolluted sinusoidal and its frequency and magnitude as must be in the 
range of IEEE and IEC standards. Nowadays, one of the main issue of PQ is harmonic pollution, this is due 
to the use of different power electronics equipment and nonlinear loads [1], [2]. A harmonic pollution 
monitoring in the distribution system is vital for the study of the root causes, harmonic level and mitigation 
[3]-[5]. This harmonic pollution can decline the PQ, it additionally can increase the power losses and lead to 
the malfunction of equipment and measurement [6], [7]. 

A comprehensive research is necessary for producing an accurate, fast and reliable method for the 
harmonic signal detection and classification [8], [9]. Numerous methods are proposed for detection and 
classification of a harmonic signal, as straightforward and fast such as Fast Fourier Transform (FFT) [10]. 
FFT is effectively employed on stationary signals [11], [12]. Yet, FFT is not an appropriate method for non- 
stationary signal analysis and contributes time information loss throughout the frequency domain conversion 
[13], [14]. To deal with non-stationary signal, short-time Fourier transform (STFT) is introduced [14], [15]. 
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Nevertheless, STFT is not proficient to detect the dynamic signal property because the window width fixed to 
affect the frequency-time resolution [6], [8], [16]. Whereas, the time-frequency resolution depends upon the 
size of the window [17], [5]. 

The constraint of STFT is resolved by utilizing Gabor transform (GT) and GT has the feature of 
extraction tool, due to the optimality concerning the time-frequency uncertainty characteristic [18]. As 
suggested by GT, the functions are well-defined as Gaussian envelope modulated by complex sinusoids with 
a fixed envelope for all frequencies [19]-[21]. However, the result shows that the method requires high 
computation complexity compared to the STFT due to the use of discrete Fourier transform (DFT) [10]. In 
order to overcome STFT and GT limitation, the wavelet transform (WT) is proposed [22]. WT can extract the 
significant information from non-stationary and it can distinguish the signal characteristics [6], [18], [23]. 
The major drawback of WT are the accuracy thoroughly relies on the chosen mother wavelet, sensitive to the 
level of noise and high computation complexity [18], [24], [25]. Subsequently, a combination of STFT and 
WT known as S-transform (ST) is presented due to mitigate the WT problem [6], [26], [27]. Furthermore, the 
ST is a reliable technique to characterize the harmonic parameters [26], [28], [22]. The ST also offers an 
extraordinary multiresolution analysis while characterizing the harmonic components [9], [18], [16], [22]. 
Based on above discussion, an evaluation and a critical review of TFDs for harmonic signal detection and 
classification are crucial and need to be done. 

The main concern of this paper is to evaluate and review the best TFDs in harmonic signal detection 
and classification technique. The assessment of the harmonic signal analysis is using the TFDs such as 
spectrogram, GT and ST are actualized in time-frequency domain [29], [30]. The performance of each TFD 
is compared and validated with regards to accuracy, computational complexity and the used memory size 
during the analysis. The performance of these TFDs are validated by detecting and classifying the 100 unique 
signals with numerous characteristics of harmonic signal accordingly to the IEEE Std. 1159-2009. The best 
TFD is chosen by determining the most accurate method, used low memory size and additionally low 
computational complexity. 


2. HARMONIC SIGNAL ANALYSIS 

There are four primary phases in the harmonic signal detection and classification analysis as shown 
in Figure 1. The Ist phase is the signal normalization and representation of time-frequency representation 
(TFR). Meanwhile, the 2nd phase is the estimation of fundamental signals parameters. Moreover, the 3rd 
phase is the signal characteristics identification and finally, the 4th phase is the harmonic signal 
classification. The harmonic signal is standardized where the voltage magnitude is changed into the per-unit 
system. Meanwhile, via TFDs, the outcome of the analysis is presented in the time-frequency domain. 

The parameters comprise of the RMS fundamental voltage, total waveform distortion, instantaneous 
of RMS voltage, total nonharmonic distortion and total harmonic distortion are calculated. Lastly, according 
to the IEEE Std. 1159-2009, this information is utilized as the rule-based classifier input due to classify the 
harmonic signals. 
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Figure 1. Flow chart of harmonic signal detection and classification 
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3. HARMONIC SIGNALS MODELING 

Due to model the signals for analysis, a few parameters for each unique signal are proposed and 
permitted to be altered according to the IEEE Std. 1159-2009. The signal model can be written as a complex 
exponential signal and well-defined as [31], 


Xua (O) =e? + Ae? (1) 


Whereby, fo is the fundamental signal frequency and f; is the harmonic or interharmonic frequency and t is 
the time, f=250 Hz, A=0.25 for harmonic and f; =275 Hz, A=0.25 for interharmonic. 


4. THE TIME-FREQUENCY DISTRIBUTIONS 

Time-frequency distributions (TFDs) are excellent methods that presented a signal in time- 
frequency representation and recognized as time-frequency representation (TFR). In the subsequent sections, 
the TFDs consist of the spectrogram, GT and ST are explained. 


4.1. Spectrogram 
The spectrogram is a method to represents a signal energy distribution in jointly time and 
frequency [32]. It is well-defined as, in this research, the Hanning window is carefully chosen as its lower 
peak side lope has a feature of narrow effect on the frequency components. For this research, the length of 
the window is 512 and the frequency and time resolution for spectrogram is fixed for all frequencies. 
2 


ee) 


P(t, f) =| | xwa Ne?" de (2) 


—0 


4.2. Gabor Transform 

The descriptor of signal’s local property of Gabor transform (GT) is obtained from a set of functions 
that are condensed in frequency and time domains [25]. The GT is well expressed by, in GT, Hanning 
window is used as well as a spectrogram, however dissimilar of window length. The resolution of frequency 
and time for GT is always same for all frequencies. 


C(n,k) = Í x(t)h' (n,k)dt (3) 


—0 


4.3. S-transform 

S-transform (ST) is a hybrid of short time Fourier transform (STFT) and wavelet transform and used 
as a time-frequency spectral localization method [17]. ST utilizes a scalable of the Gaussian window and the 
resolution of frequency based on the simultaneous localization of the real and imaginary spectra [5]. The ST 
is well written as, 


o (1-4) f? 
ST(t, f)= J rot 2 erfd (4) 
2 T 
E e20? (5) 
O T 
of)= + © 
|| 


Whereby h(t) is the signal, g(t) is the scalable Gaussian window and o(f) is a control parameter for the 
Gaussian window. ST offers superior frequency resolution for lower frequency. Meanwhile, for higher 
frequency, ST delivers good time resolution. The extraction of frequency components is achieved by utilizing 
wide window (low frequency) and narrow window (high frequency) due to compose high frequency 
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components [32]. In this research, spectrogram and GT with the utilization of Hanning window are used to 
calculate the time and frequency. Meanwhile, the scalable Gaussian window is used for ST method. 


5. SIGNAL PARAMETERS 

Harmonic signals parameters are estimated from the TFR and the signal parameters consist of 
instantaneous RMS voltage and RMS fundamental voltage, instantaneous total waveform distortion (TWD), 
instantaneous total harmonic distortion (THD) and instantaneous total interharmonic distortion (7nHD). 


5.1. Instantaneous RMS Voltage 
Root-mean square (RMS) voltage, Vins is defined as [32], 


fs 


Va = [| PE Pah D 


0 
Whereby Px(t f) is the TFR signal and fs is sampling frequency. 


5.2. Instantaneous RMS Fundamental Voltage 
From the TFR, the instantaneous RMS fundamental voltage, Virms(t), can be calculated using [31], 


Sri 


Vins(t) = [2] P.O. Pdf 
tio (8) 
Íri E fo +25 Hz 
fo = fy-25HZ 


Whereby P,(t,f) is the TFR signal, fo is the fundamental frequency. 


5.3 Instantaneous Total Waveform Distortion 
The total waveform distortion, TWD is well-defined as the relative signal energy existing at a non- 
fundamental frequency and expressed as [31], 


Viet Miva): 


TWD( jad m V O (9) 
Irms 


5.4. Instantaneous Total Harmonic Distortion 
Total harmonic distortion, THD, is utilized due to measure the harmonic content in a waveform and 
formulated as [31], 


H 


h=2 V, y rmt) i 
V, salt) 


THD(t) = (10) 


5.5. Instantaneous Total Nonharmonic Distortion 
A signal also comprises interharmonic components and the interharmonic level can be calculated 
using [31], 


Wan (t) 4 = ee Vis (t) f 
Virais (t) 


TnHD(t) = (11) 
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6. SIGNAL CHARACTERISTIC 

The characterization of signals is obtained from the calculated signal parameters. Average of total 
harmonic distortion, THDave and total nonharmonic distortion, TnHDave can be calculated from instantaneous 
total harmonic distortion, THD(t) and instantaneous total nonharmonic distortion, TnHD(t), respectively. 
These characteristics can be expressed as [31], 


1 T 
Vae =— | Vins (tat 12 
rms,ave T J ak ) ( ) 
1 T 
THD,,, = z | THD at (13) 
0 
1 T 
THD e == [TaHD (oat (14) 
0 


7. SIGNAL CLASSIFICATION 

A deterministic classification method which is the rule-based classifier extensively employed in the 
practical application. This method is easy to be implemented and its performance is much reliant on the 
threshold settings and expert rules. The flow chart in Figure 2 describes a rule-based classifier for the 
harmonic signals. Furthermore, the threshold settings are set according to IEEE Std. 1159-2009. 
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Figure 2. The rule-based classifier flow chart for harmonic signals 


8. PERFORMANCE MEASUREMENTS OF TIME-FREQUENCY DISTRIBUTIONS 
The performance and viability of proposed method depend on the accuracy of this technique and the 
evaluation of accuracy can be done as follows. 
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8.1. The Analysis Accuracy 

The analysis accuracy is calculated through the signal characteristics measurement accuracy. 
To measure the measurement accuracy, the mean absolute percentage error (MAPE) is employed as an 
accuracy index. The lower the MAPE, the better the performance of the signal characteristics 
measurement [32]. It can be written as, 


1 Alx m -x, (n) 
MAPE = i m |100% 
NÈ 0) x 100 % (15) 


Whereby x(n) is an actual value, x,,(n) is measured value and N is the data number. The smaller value of the 
MAPE, the better the accuracy is. 


8.2. The Computational Complexity of the Analysis 

The computational complexity of TFDs to represent the TFR of signals reliant on the parameters 
setting of each TFD. The spectrogram and ST both used fast Fourier transform (FFT) that deals fast 
computation. Meanwhile, GT executes discrete Fourier transform (DFT) and has more computational 
complexity compare to FFT [18]. The computational complexity index can be calculated as follows, 


CF roire = N = N, (N, log 2 (N, )) (16) 
Ci Gabor =N-N, (N?) (17) 
Cr =N-N, (N, +N, log, (N,,)) (18) 


8.3. The Memory Size of the Analysis 
The signal length number and TFD parameters setting influence the size of memory used for 
TFR [18]. The memory is expressed as, 


N,(N-N,) 
Memory jincarteD T “NO (19) 


s 


9. RESULTS AND ANALYSIS 

This section discuss the results of harmonic signals detection and classification analysis by using 
TFDs and the best TFD is chosen based on the accuracy, computational complexity and memory size during 
the TFDs analysis. This section discuss the results of harmonic signals detection and classification analysis 
by using TFDs and the best TFD is chosen based on the accuracy, computational complexity and memory 
size during the TFDs analysis. 


9.1. Harmonic Signal Analysis by Using Spectrogram 

Figure 3(a) and 3(b) depict harmonic signal in the time domain and the TFRs using spectrogram. 
The TFRs show the signal parameters comprise of the fundamental frequency at 50Hz and the 7th harmonic 
component at 350Hz, respectively. Furthermore, Figure 3(c) shows that the harmonic voltage has contributed 
to the rise of the RMS voltage from 1.0 to 1.17 pu. However, the RMS fundamental voltage still remains at 
1.0 pu. Besides that, the TWD and THD of the signals are 60% and zero TnHD are calculated and presented 
well in Figure 3(d). Hence, from the analysis show that there is no interharmonic signal is exist in the signal. 

The existence of interharmonic signal is analyzed with spectrogram and its TFR is demonstrates in 
Figure 4(a) and 4(b), separately. As demonstrates in Figure 4(b), the signal consists of fundamental frequency 
at 50 Hz and the interharmonic frequency at 375 Hz, respectively. Figure 4(c) shows that the existence of 
interharmonic voltage has increased the RMS voltage from 1.0 to 1.17 pu. The analysis also results the TWD 
and TnHD are 60% and as depicts in Figure 4(d), the THD is zero percent. Therefore, the results of the 
analysis show only fundamental and interharmonic signal exist in the signal. 
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Figure 3. (a) Harmonic signal in time domain, (b) Harmonic signal in TFR, (c) Instantaneous RMS voltage 
and RMS fundamental voltage, (d) Instantaneous total harmonic distortion, total nonharmonic distortion 
and total waveform distortion. 
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Figure 4. (a) Interharmonic signal in time domain, (b) Harmonic signal in TFR, (c) Instantaneous RMS 
voltage and RMS fundamental voltage, (d) Instantaneous total harmonic distortion, total nonharmonic 
distortion and total waveform distortion. 
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Figure 4. (a) Interharmonic signal in time domain, (b) Harmonic signal in TFR, (c) Instantaneous RMS 
voltage and RMS fundamental voltage, (d) Instantaneous total harmonic distortion, total nonharmonic 
distortion and total waveform distortion. 


9.2. Signal Analysis Using Gabor Transform 

Figure 5(a) and 5(b) present harmonic signal in time domain and the TFR by using GT. From the 
TFR, is show that the signal contains fundamental frequency at 50Hz and the 7th harmonic component at 
350 Hz. In the meantime, Figure 5(c) shows that the existence of harmonic voltage has increased the RMS 
voltage from 1.0 to 1.17 pu. The analysis outcome in Figure 5(d) shows that the TWD and THD of the signal 
at a value 60% and the TnHD is zero percent as figured in Figure 5(d). Thus, it is clearly shown that no 
interharmonic component exists in the signal. 

Figure 6(a) and 6(b) show the analysis outcome of the signal with interharmonic component and the 
TFR using GT. From the TFR, it shows that only the fundamental frequency at 50 Hz and the interharmonic 
component at 375 Hz are exist in the signal. Meanwhile, as shown in Figure 6(c), the existence of 
interharmonic voltage has increased the RMS voltage from is 1.0 to 1.17 pu. In addition, from Figure 6(d), 
the analysis also presents the value of TWD and TnHD is 60% and zero percent for THD. Thus, from the 
analysis clearly show that only fundamental and interharmonic components exist in the signal. 
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Figure 5. (a) Harmonic signal in time domain, (b) Harmonic signal in TFR, (c) Instantaneous RMS voltage 
and RMS fundamental voltage, (d) Instantaneous total harmonic distortion, total nonharmonic distortion 
and total waveform distortion 
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Figure 5. (a) Harmonic signal in time domain, (b) Harmonic signal in TFR, (c) Instantaneous RMS voltage 
and RMS fundamental voltage, (d) Instantaneous total harmonic distortion, total nonharmonic distortion 
and total waveform distortion 
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Figure 6. (a) Interharmonic signal in time domain, (b) Harmonic signal in TFR, (c) Instantaneous RMS 
voltage and RMS fundamental voltage, (d) Instantaneous total harmonic distortion, total nonharmonic 
distortion and total waveform distortion 
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9.3. Signal Analysis Using S-transform 

Figure 7(a) and 7(b) present harmonic signal in time domain and the TFR by using ST. From the 
TER, it is clearly figured that the signal comprises the fundamental frequency at 50Hz and the 7th harmonic 
component at 350 Hz. In the meantime, Figure 4(c) shows that the existence of harmonic component has 
increased the RMS voltage from 1.0 to 1.1 pu. In addition, Figure 4(d) demonstrates the magnitude of TWD 
and THD is 10% and zero percent for the TnHD. Hence, from the analysis, it is plainly shown that no 
interharmonic component exists in the signal. 
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Figure 7. (a) Harmonic signal in time domain, (b) Harmonic signal in TFR, (c) Instantaneous RMS voltage 
and RMS fundamental voltage, (d) Instantaneous total harmonic distortion, total nonharmonic distortion and 
total waveform distortion 


Figure 8(a) and 8(b) present the interharmonic signal and the TFR by using ST, respectively. The 
signal has a fundamental component at 50 Hz and interharmonic component at 375 Hz. Figure 8(c) shows 
that the interharmonic voltage contributes to the rise of RMS from 1.0 to 1.1 pu. In addition, the signal 
analysis presents the magnitude of TWD and TnHD with a value of 10% and zero percent for the THD as 
shown in Figure 8(d). Thus, the analysis clearly shows that only the fundamental and interharmonic 
components exist in the signal. 


Int J Elec & Comp Eng, Vol. 8, No. 6, December 2018 : 4603 - 4618 


Int J Elec & Comp Eng ISSN: 2088-8708 O 4613 
Power Quality Signal Time-Frequency Representation 
1.5 500 0.5 
450 0.45 
k 400 0.4 
_ 05 = od 0.35 
a = 300 
= > 0.3 
So © 250 
g A ba 0.25 
-0.5 = 0.2 
150 
i 100 0.15 
so ee 01 
-1.5 0 0.05 
20 40 60 80 100 120 140 160 180 200 220 240 20 40 60 80 100 120 140 160 180 200 220 240 
Time (msec) Time (msec) 
(a) (b) 
Q Instantaneous Total Harmonic Distortion 
Instantaneous RMS Voltage S 400 
15 o 
2 T 59 
a S 
o 9 0 
U oO 
31 E 2 40 60 80 100 120 140 160 180 200 220 240 
3 Time (msec) 
< Q Instantaneous Total Nonharmonic Distortion 
05 g” 
“20 40 60 80 100 120 140 160 180 200 220 240 Q 
Time (msec) T 50 
Instantaneous RMS Fundamental Voltage 8 © — 
1s e 0 4 60 8 100 120 140 160 180 200 220 240 
2 Time (msec) 
o Q Instantaneous Total Waveform Distortion 
| ~ 100 
z 8 
£ EE) 
< G 
05 5 
a 


0 4 
Time (msec) 


(c) 


Figure 8. (a) Interharmonic signal in time domain, (b) Harmonic signal in TFR, (c) Instantaneous RMS 
voltage and RMS fundamental voltage, (d) Instantaneous total harmonic distortion, total nonharmonic 
distortion and total waveform distortion 


9.4. The Performance Analysis of the Time-frequency Distributions 

To identify the best method for harmonic signal analysis, TFDs are compared in terms of the 
accuracy, used memory size and computational complexity of the algorithm. 100 unique signals are 
generated and employed in the analysis. 


9.4.1. The Accuracy of the Analysis 

Harmonic signals are classified and the mean absolute percentage error (MAPE) of the signal 
characteristics is calculated as shown in Table 1. The table indicates that the ST gives an excellent accuracy 
of Vins, THD and TnHD. This is due to the method offers an excellent time resolution for high frequency and 
excellent frequency resolution for low frequency. 


Table 1. MAPE of Simulation Results of the TFDs 


8 TFDs 
ae 
ge 
Bay S Spectrogram Gabor Transform S-transform 
S 
U 
VI MSave 0.1572 0.5293 0.0426 
THDave 0.1551 0.9967 0.0541 
TnHDave 0.1595 0.9331 0.0533 
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The tabulation of Table 1 is presented in the bar graph in order to have a good observation on 
MAPE value. The lower the MAPE value, the accurate the results are. Hence, the result plainly shows that 
the ST is the best method to measure the signal characteristics. 
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Figure 9. Average of MAPE for the TFDs 


9.4.2. The Computational Complexity of the Analysis 

The computational complexity of the analysis for each TFD is tabulated in Table 2. As demonstrates 
in the table, spectrogram, GT and ST contribute same computational ratio. This is because the windows 
length and a number of the signal are same for each TFD. 


Table 2. Computation Complexity of the TFDs 


' TFDs 
Signal Spectrogram Gabor Transform S-transform 
Normal 20,680,704 1,041,408,000 22,978,560 
Harmonic 20,680,704 1,041,408,000 22,978,560 


Interharmonic 20,680,704 1,041,408,000 22,978,560 


The bar graph as presents in Figure 10, point out that spectrogram offers the lowest computational 
complexity, while GT is the highest. Hence, it is summarize that the spectrogram offers the fastest speed for 
harmonic signal analysis and then ST while GT is the slowest. 


Computation Complexity 
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Figure 10. An average of computation complexity used for the TFDs 
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9.4.3. The Memory Size of Data Analysis 

Table 3 presents the average memory size used during the signal analysis for each TFD. Meanwhile, 
the graph as presented in Figure 11 show that the GT offers the lowest used memory size while spectrogram 
and ST utilized the biggest during the analysis. The memory size for spectrogram and ST are similar due to 
the same window size in the analysis. 


Table 3. Used Memory Size of TFRs for the TFDs 
TFDs, (Mbyte) 


Signal Spectrogram Gabor Transform S-transform 
Normal 2,297,856 2,250,000 2,297,856 
Harmonic 2,297,856 2,250,000 2,297,856 
Interharmonic 2,297,856 2,250,000 2,297,856 
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Figure 11. Average of memory size used for the harmonic signal analysis 


9.4.3. The Best Time-frequency Distribution Technique 

The results of the accuracy, computational complexity and used memory size of are depict in 
Figure 12. The performance evaluation of the best TFD based on the highest priority of criteria which is an 
accuracy, followed by the computational complexity and memory size. 
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Figure 12. The performance comparison of TFDs 
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As demonstrate in Figure 12, the spectrogram offers the lowest computational complexity. 
However, it used high memory size and provides the lowest accuracy. The poorest TFD is GT that request 
high computational complexity and offers low accuracy. Aimed at ST, it offers the highest accuracy, low 
memory size and additionally low computational complexity. Thus, the outcome, unmistakably, demonstrates 
that the ST is the best TFD for harmonic signal detection. 


9.4.3. The Classification of Harmonic Signal 

ST has been identified as the best TFD in this research for harmonic signals detection. The results of 
the signals classification using the ST are present in Table 4. The 100 unique signals with numerous 
characteristics were generated and classified. The classification results show that the ST provides 100% 
correct classification for all harmonic signals. From the outcomes, it can be inferred that ST is the best 
technique for harmonic signal detection and classification. 


Table 4. Performance of Power Quality Signal Classification for S-transform 
S-transform 


Signal Number of data sets % Correct Classification 
Harmonic 100 100 
Interharmonic 100 100 
Normal 100 100 


10. CONCLUSION 

The main concern of this paper is to evaluate and review the best TFD in harmonic signal detection 
and classification. The evaluation of the harmonic signal analysis using the TFDs actualized in time- 
frequency domain and compared in terms of accuracy, computational complexity and used memory size 
during the analysis. The performance of these methods are confirmed by detecting and classifying 100 unique 
signals with numerous characteristics of the harmonic signal. The outcomes demonstrate that ST is the best 
TFD with a criteria of the most accurate method, used low memory size and additionally low computational 
complexity. 
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